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1  .  OVERVIEW 


1.1  Introduction 

This  Phase  I  Technical  Report  documents  the  successful  establishment  of  the  feasibility  of  a 
novel  diaphragm  compressor  configuration  for  cryocoolers.  We  established  the  feasibility  by 
successfully  building  a  Proof  of  Principle  (PoP)  model  of  this  novel  compressor.  This  report 
presents  the  analytical  basis  and  the  results  of  an  experimental  test  program  conducted  to 
confirm  the  potential  benefits  of  the  diaphragm  compressor  design.  We  conclude  that  novel 
compressors  of  essentially  infinite  life  can  be  produced  by  this  innovative  design  concept. 


1.2  Existing  Compressor  Approaches 

A  typical  cryocooler  consists  of  a  compressor,  displacer  and  pressurized  gaseous  helium  (GHe). 
The  compressor  generates  a  pressure  wave  in  GHe  which  is  used  to  reciprocate  the  displacer  to 
produce  cold.  Its  function  is  to  efficiently  generate  the  pressure  wave. 

Space  quality  cryocoolers  should  have  extremely  long  life  and  high  reliability  to  achieve  the  space 
mission  objectives.  Ideally  long  life  requires  elimination  of  all  rubbing  parts  that  produce  friction 
and  wear.  Existing  Stirling  cryocoolers  avoid  rubbing  by  using  flexures  to  suspend  the 
reciprocating  piston.  These  pistons  use  a  very  small  clearance  to  reduce  gas  leakage.  Such 
clearance  seals  however  are  sensitive  to  manufacturing  inaccuracies.  Elimination  of  pistons, 
flexures  and  clearance  seals  will  obviously  result  in  a  compressor  that  is  free  of  wear-prone 
parts.  The  objective  of  our  diaphragm  compressor  concept  is  to  develop  such  a  design 
methodology. 
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1.3  Diaphragm  Compressor  Design 

In  our  approach,  a  diaphragm  is  vibrated  by  magnetic  fields  to  produce  the  pressure  wave.  A 
diaphragm  is  defined  as  a  thin  circular  plate  (with  no  central  hole  and  fixed  at  the  edges)  whose 
deflection  is  of  the  order  of  its  thickness.  The  magnetic  fields  are  produced  by  an  attraction 
electromagnet.  Unlike  the  existing  Oxford  design,  the  diaphragm  compressor  does  not  use 
pistons,  clearance  seals,  flexures  or  other  reciprocating  members.  The  diaphragm  sweeps  and 
seals  the  working  volume  of  the  cryocooler.  Because  this  design  eliminates  many  of  the 
precision  parts  that  are  commonly  used  in  existing  cryocoolers,  it  provides  an  added  benefit  of 
reduced  cost  and  improved  performance. 

Building  the  PoP  compressor  model  has  given  us  valuable  insight  which  makes  us  believe  that 
our  diaphragm  compressor  will  be  more  reliable.  Use  of  a  diaphragm,  however,  carries  the 
recognized  risk  of  high  cycle  fatigue  failure.  In  a  five  year  cryocooler,  the  diaphragm  accumulates 
about  1010  cycles  of  reversing  stresses.  A  diaphragm  that  can  sustain  such  a  large  number  of 
cycles  without  failure  would  require  that  these  stresses  be  low  in  magnitude.  Our  aim  is  to 
accomplish  this  by  using  large  diameter  and  small  stroke  geometries. 

This  report  documents  in  detail  the  results  of  analytical  and  experimental  programs  executed  to 
develop  the  diaphragm  compressor.  The  report  concludes  positively  that  it  is  feasible  to  design 
this  novel  diaphragm  compressor  and  recommends  moving  forward  to  Phase  II  for  developing  a 
fully  functional  engineering  prototype  compressor. 


1.4  Advantages  of  Diaphragm  Compressors 

The  design  concept  we  propose  can  be  applied  to  cryocoolers  of  the  following  types:  Stirling 
cycle,  pulse-tube,  acoustic  or  magnetic.  (It  can  also  be  applied  to  several  other  applications 
such  as  small  industrial  pumps,  refrigerator  compressors,  automotive  pumps  etc. ).  Its  primary 
function  is  to  generate  a  pressure  wave  in  a  charged  fluid  using  large  bore  and  short-stroke 
diaphragm  geometries  . 


This  diaphragm  compressor  does  not  use  the  classical  piston-in-cylinder  approach.  Hence  it 
does  not  employ  any  reciprocating  parts  such  as  pistons  or  linear  motors.  Thus  it  does  not  have 
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any  rubbing  parts.  The  elimination  of  rubbing  parts  increases  the  life  of  the  compressor. 
Diaphragm  compressors  also  offer  higher  power  densities  because  the  attraction  electromagnets 
generate  larger  force  in  a  given  volume  than  conventional  linear  motors  [1],  Higher  power 
densities  translate  to  smaller  weight  and  a  lighter  compressor. 

Diaphragm  compressors  need  fewer  precision  parts  than  the  piston  compressors.  Precise 
alignment  of  a  piston  and  a  cylinder  is  not  necessary  for  these  compressors.  Hence  they  will  be 
more  rugged  and  less  expensive. 

1.5  Phase  I  Technical  Objectives 

The  primary  objective  of  the  Phase  I  effort  is  to  establish  the  feasibility  of  this  diaphragm 
compressor.  This  effort  therefore  answers  the  critical  feasibility  questions,  namely,  will  this 
compressor  work  in  the  way  it  is  supposed  to  work,  and  can  the  diaphragm  be  sized  for  a 
space-type  cryocooler  while  maintaining  low  stress?  The  specific  technical  objectives  of  Phase  I 
are  therefore: 

o  Design  and  fabricate  a  proof  of  principle  model  of  the  compressor 
o  Determine  maximum  stresses,  and  swept  volumes  produced 

1.6  Phase  I  Technical  Approach 

The  Phase  I  began  by  a  kick-off  meeting  with  the  government  technical  personnel  wherein  the 
performance  goals  were  discussed.  During  this  meeting  it  was  agreed  that  the  construction  of  a 
proof  of  principle  model  of  the  compressor  is  the  key  requirement  to  establish  feasibility.  Based 
on  these  discussions,  a  preliminary  design  of  the  PoP  compressor  was  developed.  The  base 
line  requirements  such  as  frequency,  charge  pressures,  pressure  ratios  and  the  swept  volumes 
were  established  first  with  the  help  of  the  Project  Monitor.  These  parameters  were  then 
employed  to  analytically  evaluate  the  maximum  stresses,  deflections  and  volumes  under  the 
combined  action  of  the  attraction  and  gas  loads.  A  proof  of  principle  compressor  was  then 
fabricated  and  tested  up  to  300  Hz  to  establish  the  soundness  of  this  approach. 


The  work  was  performed  as  five  tasks  that  are  designed  to  meet  specific  technical  objectives. 
These  tasks  are: 

1 .  Determine  the  requirements  for  the  compressor 

2.  Determine  the  requirements  for  the  diaphragm  and  size  it 

3.  Analyze  the  stresses  generated  n  the  diaphragm 

4.  Design  and  fabricate  a  Proof  of  Principle  compressor 

5.  Test  the  Proof  of  Principle  cor  pressor 

The  results  of  these  five  tasks,  which  encompass  both  the  analytical  and  experimental  aspects, 
provide  a  strong  indication  that  a  diaphragm  compressor  that  can  operate  for  five  years  can  be 
developed  and  fabricated. 


(a)  Compresso:  Design  Software 

The  purpose  of  this  software  is  to  predict  the  performance  of  the  compressor  once  its  geometric 
dimensions  and  tt  e  operating  parameters  are  known.  This  will  involve  translating  the  design 
requirements  intr  a  sizing  algorithm.  The  basic  approach  is  to  develop  a  theoretical  analysis  for 
sizing  individua1  components  first,  developing  subroutines  to  code  them  and  combine  all  these 
subroutines  into  a  system  flow  chart.  This  approach  breaks  down  the  complex  task  of  developing 
a  large  unmanageable  code  into  a  manageable  number  of  smalfer  subroutines  which  can  be 
easily  coded  and  debugged. 

We  have  r  eveloped  a  baseline  compressor  design  software  as  part  of  Phase  I  effort.  We  will  be 
extendin'-  and  refining  this  software  into  a  comprehensive  stand-alone  software  module  in  Phase 
II.  Figure  T shows  a  typical  input/output  of  this  baseline  software.  The  software  inputs  the 
design  specifications  of  the  compressor  and  outputs  its  sizing  and  performance  parameters. 
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E  =  youngs  modulus,  psi 

anu=  poissons  ratio 

D=  diaphragm  diameter,  in. 

t  **  diaphragm  thickness,  in 

rO  =  central  pole's  radius,  in 

pawdc  attraction  pressure  in  work  chamber,  dc  comp.,psia 

pawac  attraction  pressure  in  work  chamber,  ac  comp.,  psi 

pabdc  attraction  pressure  in  back  chamber,  dc  comp.,  psia 

pabac  attraction  pressure  in  back  chamber,  ac  comp. ,  psi 
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WORK  CHAMBER  IS  CONTOURED 

delv  = 

.2839 

Figure  1.  Output  of  the  compressor  design  software  developed  under  Phase  I 


Sn  order  to  establish  the  Feasibility  of  the  diaphragm  compressor,  we  built  two  prototype  units. 
The  first  unit  -  shown  in  Figure  2  -  used  Permanent  Magnet  bias  to  attract  a  titanium  alloy  disk  to 
which  a  soft-steel  disk  (made  of  iron  powder  embedded  in  a  silicone  matrix)  is  attached.  This 
approach  showed  significant  promise,  but  the  diaphragm  used  was  too  stiff.  To  overcome  this 
hurdle  we  revised  our  approach  to  an  all-electromagnet  approach  in  the  second  prototype.  In 
this  prototype  -  shown  in  Figure  3  -  we  used  bias  currents  in  the  electromagnets  to  carefully 
control  the  position  of  the  diaphragm.  This  approach  proved  to  be  successful  and  has  been 
documented  in  this  report. 

Permanent  Magnet  Biased  Prototype  Compressor 

This  Prototype  used  two  NdBFe  permanent  magnet  ring,  3  in.  OD  x  1.75  in.  ID  x  0.15  in.  thick 
embedded  in  the  steel  core  of  two  electromagnets  as  shown  in  Figure  2.  The  diaphragm  is 
made  of  3.66  in.  diameter  Titanium  alloy  disk  which  is  0.05  in.  thick  to  which  a  flexible  steel  disk, 
which  is  3.60  in.  diameter  and  0.2  in.  thick  is  attached  as  shown  in  Figure  2.  The  electromagnet 
coil  is  mounted  in  an  annular  slot  with  4.66  in.  diameter  as  shown.  The  diaphragm  is  held  in 
position  by  precisely  machined  annular  rings  made  of  non-magnetic  stain  less  steel.  The  two 
electromagnets  are  bolted  together  by  a  set  of  6  bolts. 

Figure  2  also  shows  the  flexible-steel  disk  that  is  fabricated  during  this  prototype  building  activity. 
Realizing  that  the  thickness  of  the  disk  must  be  large  enough  to  allow  closure  of  flux  path,  we 
developed  this  soft-steel  disk  as  part  of  this  Phase  I  activity.  This  flexible  disk  is  made  by  mixing 
fine  iron-powder  into  a  silicone  mould  and  using  an  appropriate  catalyst  to  embed  the  iron- 
powder  particles  in  the  silicone  matrix.  The  silicone  matrix  offered  the  flexibility  that  is  needed  by 
the  diaphragm  while  the  iron  powder  provided  the  medium  for  the  flux  path.  This  approach 
resulted  in  a  new  material  aptly  called  flexi-steel  disk  since  it  combined  the  magnetic  properties 
of  rigid  steel  and  flexibility  of  silicone. 

During  the  testing  of  this  prototype  however,  we  found  that  the  Titanium  alloy  diaphragm  is  too 
thick  and  hence  the  electromagnet  is  not  strong  enough  to  vibrate  the  diaphragm.  As  a  result, 
for  the  purpose  of  Phase  I  activity  we  elected  to  make  another  prototype  which  used  bias 
currents  to  position  the  diaphragm. 
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PARTS  LIST! 

1- COIL  ASSEMBLY 

2- PERMANENT  MAGNET  RINGS 

3- IRON  POVDER-SILCONE  DISC 

4- AIRPAX  PASSIVE  MAGNETIC  PICKUP  70085-1010- 

5- RIGHT  HOUSING 

6- LEFT  HOUSING 

7- PCB  PRESSURE  TRANSDUCER  101A06 

8- SEALING  RINGS  LEFT  AND  RIGHT 

9- PARKER  O-RINGS 

10- 10-32  X  .25  LONG  STN.  STL.  SOCKET  HEAD  CA 

11- AIR  INLET-0.125  TUBING  SILVER  SOLDERED 
TO  HOUSING.  NUPRO  CHECK  VALVE  SS-2C-1/3 

12- OUTLET  FOR  COIL  VIRES 

13- TIT AN1UM  DISC-E.B.  WELDED  TO  SEAL  RING 


Figure  2.  Permanent  Magnet  Bias  Compressor  prototype.  This  is  the  first  prototype  that 
was  built,  but  the  diaphragm  that  was  used  proved  to  be  too  stiff 
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All-electromagnet  biased  compressor  prototype 


The  second  prototype  that  we  built  used  a  single  electromagnet  to  position  the  diaphragm  and 
vibrate  it,  thus  creating  the  pressure  wave  needed.  This  prototype  is  described  in  more  detail  in 
Section  2  and  in  the  later  parts  of  this  report.  Essentially  this  prototype  used  two  electromagnets 
sandwiching  a  diaphragm  and  vibrating  it.  We  attached  two  steel  disks  to  the  diaphragm  in  order 
to  close  the  flux  circuit  and  also  to  generate  larger  forces.  This  approach  proved  to  be 
successful.  This  Phase  I  report  documents  the  design  methodology  used  to  design,  fabricate 
and  test  this  type  of  compressor. 
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Figure  3.  All  Electromagnet  Compressor  Prototype.  This  is  the  second  prototype 
compressor  that  was  built  under  this  Phase  I  effort.  This  prototype  has  been  successfully 
used  to  generate  the  desired  pressure  wave 


The  unique  features  of  this  second  prototype  are  the  use  of: 

-  a  three  point  seal  support  of  the  diaphragm  to  increase  its  life 

-  a  switching  amplifier  to  reduce  power  loss 

-  a  low  inertia  moving  member  to  reduce  force  required 
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As  shown  in  Figure  3,  this  prototype  compressor  used  a  single  diaphragm,  made  of  spring  steel, 

5  in.  diameter  and  0.020  in.  thick  which  is  supported  on  O-  rings  at  4  in.  diameter.  The 
diaphragm  used  a  disk  of  3.2  in.  OD  and  0.125  in.  The  disk  enhances  the  attraction  force 
produced  on  the  diaphragm.  The  disk  is  mounted  on  a  10-32  screw  and  spaced  by  a  plastic 
shim  which  is  0.015  thick  in  diameter.  The  thread  of  this  screw  is  scratched  in  order  to  equalize 
the  pressures  in  the  compression  and  back  chamber.  Initially  we  tried  a  three  point  support  for 
the  diaphragm  and  found  that  it  caused  bending,  thus  centering  it  was  difficult.  We  then 
replaced  it  by  the  two  point  seal  support,  with  the  seals  spaced  on  either  side  of  the  diaphragm  at 
the  same  diameter.  Our  tests  indicated  that  ensuring  symmetric  gaps  is  the  key  to  the 
successful  implementation  of  the  prototype.  The  electromagnets  used  coils  in  a  slot  which  has  2 
in.  ID  x  2.8  in.  OD  and  0.38  in.  height.  The  coils  are  made  of  gage  22  magnet  wire  with  120 
turns.  The  cylindrical  face  is  contoured  to  take  up  the  shape  of  the  deflected  diaphragm.  The 
electromagnet  is  driven  by  an  external  switching  amplifier  using  a  laboratory  signal  source.  A 
photonic  vibration  sensor  and  a  pressure  transducer  are  used  to  monitor  the  vibrations  of  the 
diaphragm  and  the  pressure  wave  produced  by  it  respectively. 

While  developing  the  prototype  we  investigated  two  options  of  electromagnet  configurations  and 
three  options  of  its  drive  circuit.  The  electromagnets  can  be  arranged  in  NSSN  configuration  so 
that  they  repel  each  other  (the  flux  from  one  electromagnet  loops  back  to  itself)  or  in  the  NSNS 
configuration  so  that  they  attract  each  other  (the  flux  from  one  electromagnet  goes  through  other 
electromagnet  before  returning  ).  We  found  that  these  configurations  options  made  no  difference 
in  pressure  output.  We  also  tried  three  different  options  of  driving  the  electromagnets  -  first  with 
no  capacitor,  second  with  a  capacitor  on  one  electromagnet  and  a  third  with  capacitor  on  both 
electromagnets.  We  found  that  the  placement  of  the  capacitor  did  make  difference  in  the 
pressure  output.  The  option  in  which  the  capacitor  is  connected  across  one  electromagnet 
yielded  the  greatest  pressure  output  while  the  other  two  options  yielded  poor  results.  Hence  we 
chose  this  drive  circuit  option. 

After  gaining  a  thorough  understanding  of  the  operating  principle  of  this  prototype,  we  started 
investigating  the  maximum  pressure  output  that  the  prototype  can  deliver.  Our  initial  tests 
however  indicated  that  the  maximum  pressure  output  strongly  depends  on  the  dead  space  in  the 
compression  chamber.  During  this  prototype  building  activity  we  did  not  have  full  control  on  the 
dead  space,  and  hence  the  performance  is  lower  than  what  can  be  expected.  Even  with  this 
deficiency,  we  have  successfully  vibrated  the  diaphragm  at  up  to  300  Hz  and  developed  pressure 
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amplitudes  of  up  to  3  psi  at  atmospheric  pressure.  We  then  filled  the  compressor  with  air 
pressurized  at  20  psi;  in  this  charged  mode,  we  could  vibrate  the  diaphragm  at  frequencies  of 
up  to  60  Hz  and  generate  pressure  waves  of  up  to  5  psi  in  amplitude.  Thus  we  showed  that 
even  this  preliminary  prototype  can  yield  pressure  ratios  of  up  to  1.7. 

In  summary,  we  achieved  pressure  ratios  of  up  to  1 .7  at  frequencies  of  up  to  60  Hz.  Since  the 
space  quality  cryocoolers  require  pressure  ratios  of  2  and  frequencies  of  40  Hz,  the 
demonstrated  parameters  are  within  the  realm  of  the  cryocoolers.  Hence  we  believe  that  this 
technology  has  high  potential  and  can  be  refined  in  Phase  II  and  integrated  into  the  cryocooler. 
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2.  DIAPHRAGM  COMPRESSOR  CONCEPT 

2.1  Compressor  Overview 

Figure  4  shows  a  schematic  of  this  compressor.  The  compressor's  primary  function  is  to 
generate  a  pressure  wave.  It  uses  large  bore  and  short-stroke  diaphragm  geometries.  It 
comprises  of  the  mechanical  compressor  itself  and  its  electronic  controller.  The  mechanical 
compressor  is  connected  to  the  payload  volume  of  a  cryocooler.  The  compressor  hardware  is 
driven  by  the  electronic  controller.  The  electronic  controller  is  a  separate  unit  that  is  attached  by 
cables  to  the  mechanical  compressor.  It  houses  the  signal  generator  and  power  amplifiers  which 
produce  the  oscillatory  current  needed  for  the  compressor. 

The  compressor  operates  as  follows.  Two  attraction  electromagnets  vibrate  the  diaphragm  in 
pull-pull  mode,  thus  generating  the  required  pressure  wave.  This  pressure  wave  is  transmitted 
through  the  charged  fluid  into  the  displacer  volume  of  the  cryocooler. 

Figure  5  shows  the  mechanical  compressor;  it  comprises  of  a  vibrating  diaphragm,  compression 
chamber  and  the  cylinder.  The  diaphragm  is  firmly  attached  to  the  cylinder  head  at  the  outer 
periphery.  The  cylinder  head  is  contoured  and  contains  electromagnets  (comprising  a  coil  and  a 
core).  The  diaphragm  is  separated  from  the  cylinder  head  by  a  small  compression  volume.  A 
small  gap  between  the  diaphragm  and  the  electromagnet  thus  constitutes  a  compression 
chamber  in  which  the  diaphragm  vibrates.  This  volume  contains  the  pressurized  gas  that  is 
alternately  compressed  and  expanded  by  the  diaphragm.  These  compressions  and  expansions 
create  the  desired  pressure  wave. 

The  electronic  controller  comprises  of  a  signal  generator  and  a  power  amplifier.  It  is  energized 
by  a  DC  power  supply,  if  necessary,  a  closed  loop  controller  can  be  added  to  accurately  control 
the  position  of  the  diaphragm  with  respect  to  the  current  signal.  The  signal  generator  and  power 
amplifier  are  chosen  so  that  their  frequency  and  power  capacity  are  commensurate  with  the 
application  requirements. 
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Figure  6  shows  how  the  key  components  of  the  compressor  work  together.  The  signal  generator 
produces  a  sinusoidal  voltage  signal  at  a  desired  frequency.  The  power  amplifier  uses  this  signal 
to  produce  a  sinusoidal  current.  This  current  is  fed  into  the  coils  of  the  electromagnet,  thus 
generating  ampere-turns.  The  soft-iron  core  converts  the  ampere-tums  into  flux  density.  The  flux 
density  in  the  gap  produces  an  attraction  force.  This  force  causes  the  diaphragm  to  vibrate. 
These  vibrations  create  the  swept  volume.  The  swept  volume  in  the  confined  spaces  of  the 
compression  chamber  results  in  the  desired  pressure  wave.  We  give  below  the  salient  features 
of  the  key  components. 


Cylinder.  The  cylinder  comprises  of  a  contoured  head  and  attraction  electromagnets.  It  provides 
the  working  volume  needed  for  the  compressor.  The  contoured  head  permits  expulsion  of  all  the 
fluid  from  the  compression  chamber  during  the  compression  stroke  of  the  diaphragm.  The 
contoured  cavity  is  sealed  by  the  diaphragm  to  prevent  any  leakage  paths.  A  resin  mold  is 
applied  to  the  coil  to  prevent  gas  seepage.  The  dead  volume  can  be  minimized  by  contouring 
the  cavity  using  NC  machining  to  achieve  a  large  pressure-wave. 
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Figure  5.  Electromagnetic  Diaphragm  Compressor  Configuration 

Diaphragm.  The  diaphragm  is  the  critical  component  that  sweeps  the  working  volume  to 
generate  the  pressure  wave.  It  is  made  of  high-fatigue  strength  material  such  as  spring  steel  [3] 
since  its  primary  function  is  to  demonstrate  that  the  compressor  can  operate  with  long  life.  The 
diaphragm's  surface  is  finished  smooth  and  is  free  of  nicks  and  scratches  to  reduce  surface 
stresses.  The  thickness  variations  are  meticulously  controlled.  It  is  equi-spaced  between  the 
compression  and  back  chambers.  At  the  outer  edge,  the  diaphragm  is  fixed  to  the  cylinder  head. 
A  special  fixturing  design  prevents  fretting  failure.  This  design  uses  two  to  four  O-rings.  The  O- 
rings  are  pre-loaded  to  seal  the  diaphragm.  They  are  housed  in  a  circular  groove  so  that  as  the 
diaphragm  vibrates,  no  part  of  the  diaphragm  touches  the  housing.  This  special  design  feature 
prevents  fretting  failure. 

Electromagnet.  The  electromagnets  generate  the  pull-force  on  the  diaphragm.  Each 
electromagnet  comprises  of  the  coil  and  the  core.  The  coil  is  of  drop-in  design  and  does  not 
contain  a  bobbin  to  improve  the  ampere-turns.  It  is  vacuum-impregnated  with  high  temperature 
epoxy  and  is  insulated  from  the  core.  It  is  hi-pot  tested  for  insulation  breakdown  and  then  potted 
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to  prevent  any  gas  leakage.  To  seal  it,  a  high  quality  epoxy  is  used.  The  excess  epoxy  is 
removed  by  machining.  The  magnet  wire  leads  are  protected  against  rough  handling.  The  core 
is  made  of  low-carbon  steel  and  is  over-sized  to  prevent  saturation.  The  electromagnet  has  a 
central  pole  and  an  outer  pole.  Our  investigations  found  that  the  outer  pole  is  not  effective  in 
displacing  the  diaphragm  as  it  is  too  close  to  the  fixed  edge  of  the  diaphragm,  As  such  its  effect 
will  be  neglected  in  the  design  of  the  compressor. 

Controller.  The  controller  provides  the  drive  current  needed  by  the  electromagnet  coils.  A  BK 
Precision  301  IB  signal  generator  was  used  as  the  signal  source  in  Phase  I.  This  signal 
generator  can  be  replaced  by  alternative  commercially  available  signal-generating  chips  -  this 
approach  will  be  investigated  in  Phase  II  for  optimal  performance.  This  generator  can  supply  the 
required  sine  signal  at  frequencies  of  up  to  1000  Hz  or  more.  The  power  amplifier  is  made  by 
Advanced  Motion  Controls,  model  3QA20-AC.;  it  is  pulse  width  modulated, 


DEAD  VOLUME 


Figure  Block  Diagram  Showing  Interrelation  of  Key  Components . 


with  33  kHz  switching  frequency.  It  has  peak  power  output  of  5.1  kW  (  +/- 170  V  at  30  A), 
bandwidth  of  2.5  kHz  and  gain  of  8  AN  ,  and  weighs  2.5  lb.  This  current  amplifier  outputs  the 
current  signal  from  a  given  voltage  signal.  In  the  second  prototype  this  amplifier  will  be  replaced 
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by  an  optimized  unit  for  minimal  power  consumption.  Its  5.1  kW  output  power  capacity  is  ample 
and  can  easily  drive  the  30  watts  coils.  Its  gain  of  8  amp/volt  is  large  enough  to  apply  the  desired 
sinusoidal  amperes  into  the  coil.  Its  frequency  bandwidth  of  DC  to  2500  Hz  spans  the  near-DC 
frequencies  needed  for  the  diaphragm. 
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The  objective  of  the  Phase  i  program  is  to  show  that  the  diaphragm  compressor  is  realistic  and  is 
of  value  to  the  space  cryocooler  community.  In  our  design  a  circular  diaphragm,  vibrated  by  an 
electromagnet  is  used  to  produce  meaningful  pressure  waves  while  sustaining  low  level 
stresses.  To  provide  a  design  basis,  we  first  developed  a  design  methodology  which  is 
described  in  this  section. 

In  Phase  I  we  selected  a  1  W/80  K  cryocooler  as  the  target.  This  specification  encompasses  the 
space  quality  cryocoolers  as  well  as  tactical  cryocoolers.  The  latter  type  are  required  in  large 
numbers  and  current  designs  exhibit  significant  problems  with  reliability,  life  and  cost.  We 
therefore  believe  that  we  can  cover  the  space  applications  as  well  as  commercial  applications  by 
this  specification.  Even  though  the  capacity  chosen  was  relatively  arbitrary,  our  analysis 
indicates  that  we  can  successfully  scale  it  up  over  a  range  of  cooling  capacities. 


The  design  of  the  compressor  begins  by  specifying  its  performance  parameters.  These 
parameters  are  :  (1)  the  output  PV  power  Ppy,  (2)  the  efficiency  r)  ,  defined  as  the  ratio  of  the 
output  PV  power  and  the  input  electrical  power  Pln  (3)  the  operating  frequency  f  (4)  the  mean  fill 
pressure  pm  (5)  the  pressure  ratio  r,  defined  as  the  ratio  of  maximum  pressure  to  the  minimum 
pressure  and  (6)  the  phase  angle  $  by  which  the  pressure  wave  leads  the  diaphragm's  central 
deflection. 

To  generate  the  pressure  wave,  the  compressor  consumes  electrical  power  and  outputs 
mechanical  PV  power.  The  ratio  of  the  output  PV  power  and  the  electrical  power  is  called  the 
efficiency.  Obviously,  a  compressor  should  be  as  efficient  as  possible.  Existing  piston-in¬ 
cylinder  designs  usually  are  70  %  efficient.  Thus  to  output  a  typical  20  watts  of  PV  power,  its 
coils  consume  almost  30  watts  of  electrical  power. 


The  efficiency  of  electromechanical  power  conversion  depends  on  how  well  the  compressor 
chamber  is  tuned  with  the  "load"  (i.e.,  displacer).  The  compressor  is  like  a  battery.  The  displacer 
is  the  load  which  draws  power  from  this  battery.  The  impedance  of  the  displacer  (the  load)  must 
therefore  be  matched  with  the  compressor's  (battery's)  impedance  for  optimum  power  transfer. 
Knowledge  of  the  displacer  volume  is  therefore  essential  for  an  optimal  design  of  compressor. 
However,  since  the  thrust  of  this  phase  I  effort  is  to  establish  the  operating  principle,  we  have 
not  investigated  the  effect  of  the  displacer  load  on  the  compressor  design. 

Recently  a  LOM  (Linear  Oscillating  Motor)  compressor  using  an  annular  diaphragm  i.e.,  a 
diaphragm  with  a  central  hole  was  developed  [1],  Note  that,  in  contrast,  PREMAG's  design  uses 
a  diaphragm  without  the  central  hole.  The  elimination  of  the  central  hole  avoids  stress 
concentration  and  hence  increases  life.  Further,  the  prior-art  design  used  a  Linear  Oscillating 
Motor  which  uses  a  complex  flux  circuit;  in  contrast  we  used  attraction  electromagnet  which  use 
a  simplified  flux  circuit.  To  drive  the  LOM  the  prior  art  used  two  signal  generators  and  two  power 
amplifiers;  in  contrast  we  use  only  one  signal  source  and  one  power  amplifier.  A  typical 
pressure  wave  produced  by  the  prior-art  compressor  is  shown  in  Figure  7.  The  design 
parameters  of  the  PoP  model,  derived  from  this  Figure,  are: 

pm  =Mean  pressure,  ps  i  =25  psia 

r  =  Pressure  ratio  =  2 

f  =  Frequency,  cps  =  25  Hz 

s  =  Stroke,  mm  =  0.016  in.  (0.4  mm) 
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We  may  add  that  in  Phase  II  we  propose  to  develop  a  well  engineered  compressor  that  can  drive 
a  space-quality  cryocooler.  Typical  requirements  for  such  phase  II  compressor  will  be  [for  a  0.8 
W/80  K  cryocooler,  see  ref.  [6]: 


Ppv  =Output  PV  power,  W 
H  =  Compressor  efficiency 
pm  =Mean  fill  pressure,  psi 
r  =  Pressure  ratio 

f  =  Frequency,  cps 

AV  =Swept  volume,  mm3 


=  20  watts 
=  70  % 

=  150  psi  (1  MPa) 


=  40  Hz 

=  0.12  in3  (2  cc) 


5*] 


TBD 


In  addition  to  the  above,  an  optimally  designed  compressor  requires  that: 

(a)  The  diaphragm  must  have  infinite  fatigue  life  (i.e.,  maximum  stress  <  10,000  psi). 

(b)  The  compressor  must  consume  minimal  power  (i.e.,  n  >  70  %)  . 

(c)  The  compressor  must  be  light-weight  (i.e.,  <  6.6  lb). 

(d)  The  vibrations  of  the  compressor  must  be  minimal  (<.2  N  or  0.05  lb). 

In  the  next  few  sections  we  discuss  the  implications  of  these  specifications. 


3.2  Output  PV  Power 


The  output  PV  power  equals  product  of  the  area  of  the  PV  diagram  and  the  frequency  f : 

PpvmfipdV 


where  f  =  operating  frequency,  p  =  pressure  and  v  =  volume.  The  integral,  which  represents  the 
area  of  the  PV  diagram,  can  be  estimated  from  specified  pressure  drop  2p0  and  swept  volume  Av 
using: 

J pdV  =  b  (2 P0  A  V) 


where  b  is  a  PV  diagram  factor.  Typical  values  of  this  factor  are  determined  from  two  typical  PV 
diagrams  shown  in  Figure  8.  These  diagrams  are  for:  (a)  for  a  5W/65K  cryocooler  and  (b)  a 
45W/80K  cryocooler  respectively.  From  these  we  found  that  b  varies  from  33%  to  55%  [7-8]. 

We  can  use  these  values  as  empirical  range  to  estimate  swept  volume  needed.  For  example, 
for  the  phase  II  compressor,  p0  =  50  psi,  f  =  40  Hz,  Ppy  =  20  W  and  b  =  0.33  yields  Vs  =  0.14  in3, 
which  is  close  to  0.12  in3  swept  volume  specified. 
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3.3  Efficiency 


A  large  portion  of  the  power  inputted  is  used  to  produce  the  output  PV  power.  The  efficiency  n 
is  defined  as  the  ratio  of  the  output  PV  power  to  the  input  power.  The  power  inputted  into  the 
compressor  can  be  computed  from 


inp 


p 


PV 


r\ 


For  example,  a  compressor  with  Ppy  =  20  W  and  n  =  0.7  yields  Pinp  =  29  watts.  The  difference 
between  input  power  and  output  power  (i.e.,  9  watts)  constitutes  a  loss.  This  loss  obviously 
should  be  minimized.  These  power  losses  are  broadly  divided  into  the  l2R  losses  in  the  drive 
coils,  losses  due  to  leakage  of  flux,  core  losses  and  hysteresis  losses.  Approaches  to  minimizing 
these  power  losses  will  be  investigated  in  Phase  II. 


3.4  Pressure  Wave 


The  pressure  wave  is  specified  usually  by  the  mean  pressure  pm  and  the  pressure  ratio  r.  It  is 
also  expressed  in  terms  of  other  characteristic  parameters  [Marquardt  1993]  using  following 
formula: 


Pc  - 


r- 1 
M 


Pm 


r 


=  Pm'  P 
Pm 

P  n 


-  Pm-  P 


P  ■ 

r  min 


o 

o 


where 


r  =  pressure  ratio 

p0  =  pressure  wave  amplitude 

pm  =  mean  charge  pressure 
pmax  =  maximum  pressure 
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minimum  pressure 


Pmin 


For  the  model  compressor,  pm  =  25  psia  and  r  -  1.5  yields 


p0  =5  psi 

Pmax  =  30  psi 
Pmin  =  20  psi 
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4.  DIAPHRAGM 


The  key  to  successful  design  of  a  diaphragm  is  to  insure  that  the  largest  stress  concentration  is 
smaller  than  the  smallest  allowable  fatigue  stress.  The  diaphragm  in  Phase  I  was  limited  to  a 
uniform  flat  design  of  constant  thickness.  Good  insight  into  the  effect  of  diaphragm's  material 
and  geometric  parameters  on  the  performance  of  the  compressor  is  afforded  by  this  simple 
geometry.  This  simplification  allowed  us  to  develop  stress  minimization  schemes  within  the 
Phase  I  funding.  The  use  of  complex  shaped  variable  thickness  diaphragms  could  yield  further 
reductions  in  stresses  and  improve  the  life,  but  this  will  be  accomplished  in  Phase  II  of  the 
project. 


4.1  Design  Requirements 

The  following  compressor  design  requirements  will  impact  the  design  of  a  diaphragm: 


pm  =  mean  pressure 
r  =  Pressure  ratio 
max/min./amp  pressure 
f  =  Frequency,  cps 
s  =  Stroke,  mm 
V  =  swept  volume 


=  25  psia 
=  1.5 

=  30/20/5  psi 
=  25  Hz 

=  0.016  in.  (0.4  mm) 
=  0.06  in3  (1  cc) 


In  addition,  a  dominant  requirement  is  that  the  diaphragm  should  not  fail  while  generating  the 
pressure  wave.  An  ideal  diaphragm  should  therefore  : 


o  develop  least  cyclic  stresses 

o  use  minimum  attraction  force  to  generate  largest  pressure  wave 
o  be  stiffer  than  the  negative  stiffness  of  electromagnets 
o  be  compatible  with  the  working  fluid  of  the  cryocooler 


For  the  prototype  compressor,  we  size  the  diaphragm  to  meet  following  requirements: 
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develop  a  maximum  stress  of  ±10,000  psi  in  the  diaphragm 
o  use  least  attraction  force  to  generate  a  pressure  wave  of  ±5  psi 
o  be  stiffer  than  the  negative  stiffness  of  electromagnets 
o  be  compatible  with  gaseous  helium 

These  requirements  influence  the  stresses  developed  in  the  diaphragm  and  swept  volumes 
produced.  The  design  of  diaphragm  involves  determining: 

D  =  diameter  of  the  diaphragm,  in 
t  =  thickness  of  the  diaphragm,  in. 

E  =  Young's  modulus  of  diaphragm  material,  psi 
omax  =  Maximum  fatigue  stress  of  diaphragm  material,  psi 

In  the  following  we  assume  that  the  diameter  is  given  a  priori  as  D  -  3.9  in.  The  problem  then 
reduces  to  finding  t  and  E  that  can  generate  the  specified  ±  5  psi  pressure  wave  without  fatigue 
failure. 

4.2  Load  Analysis 

The  diaphragm's  performance  is  determined  by  the  type,  magnitudes  and  phases  of  loads  and 
their  area  of  application.  These  loads  are  of  two  types  -  attraction  loads  and  gas  loads,  and 
originate  from  the  compression-chamber  and  back-chamber.  The  gas  pressure  is  always 
compressive  and  acts  over  the  entire  diaphragm  area  while  the  attraction  pressures  are  always 
tensile  and  act  only  over  the  pole  area.  In  reality  these  loads  vary  with  the  displacement  of  the 
diaphragm.  Hence  they  have  a  spring  characteristic  (force  varies  with  displacement)  instead  of 
being  constant  (force  is  independent  of  displacement).  We  analyze  this  displacement 
dependency  herein. 

We  broadly  classify  these  loads  as  DC  (or  static)  loads  and  AC  (or  dynamic)  loads.  The  DC 
loads  are  constant,  and  do  not  change  sign.  They  are  considered  positive  if  the  diaphragm 
deflects  out  of  the  compression-chamber.  The  AC  loads  are  variable  and  change  sign;  they  are 
characterized  by  an  amplitude  and  frequency.  These  are  defined  by  the  following  notation: 
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Attraction  Loads: 

-  Adcc  -  Aacc  sin  cot  =  Attraction  electromagnet  in  compression-chamber 
+Adcb  -  Aaob  sin  cot  =  Attraction  electromagnet  in  back  chamber 


Gas  Loads: 

Gdoc  +  Gacc  sin  (cot  -  $)  =  Gas  load  in  compression-chamber 
-  Gdcb  +  Gacb  sin  (cot  -  $)  =Gas  load  in  back  chamber 

where  the  symbol  A  denotes  attraction  load,  G  denotes  gas  load,  c  denotes  compression 
chamber  and  b  denotes  back  chamber. 

4.2.1  DC  Load 

The  DC  load  on  the  diaphragm  is  therefore  given  by: 

Pdc  =  -Adoc+Adcb+Gdcc-Gdcb 

The  attraction  and  gas  loads  are  estimated  from  analysis  developed  below.  Note  that,  because 
these  loads  are  distributed  over  different  areas,  the  diaphragm  may  deflect  even  if  the  attraction 
pressure  numerically  equals  gas  pressure.  In  Phase  I  we  use  the  double-acting  electromagnet 
approach  wherein  the  DC  loads  from  the  back  chamber  cancel  those  from  compression 
chamber.  This  is  accomplished  by  gas  pressure  equalization  and  bias  current  equalization.  In 
Phase  II  we  will  examine  the  merits  of  single  acting  versus  double  acting  diaphragm  approach. 


4.2.2  AC  Load 

The  AC  forces  from  the  compression  and  back  chamber  add  to  each  other.  The  ac  load  is 

pac  =  -  [Aacc  sin  cot  +  Aacb  sin  cot  -(Gac°  sin  (cot  -  cj>)  +  Gacb  sin  (cot  -  $))  ] 

The  AC  attraction  load  should  be  sized  to  overcome  not  only  the  usable  pressure  wave  in  the 
compression-chamber  but  also  that  in  the  back  chamber.  The  design  constraint  is  that  the 
smallest  AC  attraction  pull  must  be  used  to  get  the  largest  diaphragm  deflection  while  minimizing 
power  consumption  and  stresses  in  the  diaphragm. 
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(a)  Attraction  Load  Analysis 


The  mechanism  by  which  an  electromagnet  generates  the  attraction  pull  is  shown  in  Figure  9. 
The  attraction  magnet  consists  of  a  coil  of  N!  ampere  turns  wound  around  the  central  core  of  the 
electromagnet.  The  coil  generated  flux  has  a  tendency  to  close  the  loop  around  itself  using  the 
shortest  path  that  offers  smallest  resistance.  The  flux  emanates  from  the  coil  passes  through 
central  pole  into  the  air  gap  and  thence  into  the  diaphragm.  It  then  flows  radially  outward  in  the 
diaphragm  up  to  the  fixed  edge  and  thence  completes  the  flux  circuit  through  the  clamped 
portion.  The  core  and  the  coil  transform  the  electrical  energy  into  the  magnetic  energy.  The 
magnetic  energy  is  concentrated  in  the  air  gap  of  the  central  pole  gap  area.  In  this  gap  the 
magnetic  energy  is  transformed  again  into  mechanical  energy.  Because  the  diaphragm  takes  a 
cup-shape  after  deflection,  this  gap  varies  with  radial  position.  But  for  our  first  order  analysis,  we 
assume  it  to  be  uniform.  Analysis  of  a  simple  magnetic  circuit  -  ignoring  the  leakage  flux  - 
shows  that  the  flux  in  this  gap  is  given  by 

A„  N1 
4>  .  -LiL_£ — 

g  -  * 


where 

<f)  =  gap  in  the  flux  circuit 
Ap  =  area  of  the  central  pole 
Nl(t)  =  ampere  turns 

g  =  mean  gap  between  the  diaphragm  and  the  central  pole 

x(t)  =  instantaneous  deflection  of  diaphragm,  positive  if  into  compression-chamber 

The  diaphragm  is  sized  to  carry  the  flux  produced  by  the  two  attraction  magnets.  Hence 
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Figure  9.  Method  of  Generation  of  Attraction  Forces. 
the  gap  flux  density  and  diaphragm  flux  densities  are  given  by 


A„  g  -  x 
p  * 


B/dm2B/g  Bd 


2A. 


Bd  <  B  s at 


g  -  X 


where 

Ad  =  outer  peripheral  area  of  the  diaphragm 
Bd  =  flux  density  in  the  diaphragm 


The  attraction  force  is  governed  by 


Bl 

A(t)  .  -i-  a 

; 


where  =  gap  fiux  density  (<t>/Ap).  W®  suppose  that  the  current  contains  both  DC  and  AC 
components,  I  =  !dc+  lac  sincot,  so  that  attraction  force  can  be  written  as 


m  - 


2  g  1  -  */* 


When  AC  currents  are  small  relative  to  the  DC  currents  and  when  diaphragm  displacement  is 
small  relative  to  gap,  w®  can  rewrite  the  attraction  force  in  terms  of  DC  and  AC  components  as 


m  -  Adc 


A  „„ 

ac 


A 


dc 


2  g 


where 


denotes  the  DC  attraction  force.  Note  that  this  DC  force  increases  with  the  square  of  the  ratio 
DC  ampereturns/mean  gap  and  varies  directly  with  the  pole  area.  It  represents  an  upper  bound 
to  the  AC  component  of  attraction  force  which  itself  is  given  by 
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k  x  where 

X 


The  first  equation  above  indicates  that  the  AC  attraction  force  depends  not  only  on  the  AC 
currents  but  also  on  the  displacement.  The  first  term,  called  the  current  component,  is  created 
by  the  applied  currents  and  hence  requires  external  power.  Its  strength  depends  on  the  current 
stiffness  k,and  strength  of  AC  currents. 


The  second  term,  called  the  deflection  component,  is  caused  by  the  negative  stiffness  kx  (i.e., 
gradient  of  attraction  force  with  displacement).  This  force  is  characteristic  of  the  attraction 
electromagnet  and  does  not  require  external  energy.  As  the  diaphragm  moves  closer  to  the 
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centra!  pole,  the  gap  reduces  and  hence  the  diaphragm  experiences  greater  attraction  force. 
This  attraction  gradient  force  aids  the  process  of  generation  of  pressure  wave  without 
consuming  electrical  power.  This  peculiar  property  of  attraction  electromagnets  thus  reduces  the 
electrical  power  required  to  generate  a  specified  pressure  wave.  As  a  result,  the  efficiency  of 
electromagnetic  compressors  can  be  higher  than  that  of  piston-type  compressors.  A  simple 
Fortran  program  was  written  to  compute  the  attraction  load  using  this  analysis. 

(b)  Gas  Spring/load  analysis 

Initially  the  active  volume  of  compressor  and  the  displacer  is  filled  with  He  at  a  mean  pressure 
Pm.  The  net  volume  of  compressor  and  displacer  is  denoted  by  V.  Let  the  mean  pressure  in  the 
back  chamber  be  p'm  =  pm  and  its  net  volume  is  V  .  A  displacement  x  of  the  diaphragm  into  the 
compression  chamber  translates  into  a  pressure  wave  as  described  below. 

The  instantaneous  displacement  x  increases  the  pressure  in  the  compression  space  from  pm  to 
p(x)  and  reduces  its  volume  from  v  to  V(x).  The  volume  of  the  conical  bulge  shape  of  the 
deflection  curve  is  approximated  by  Ax/3,  (where  A  =  area  of  the  diaphragm)  so  that  the 
instantaneous  volume  reduces  to 

V(x )  -  AL  -  Ax/l  -  AL(1  -  xBL ) 

The  adiabatic  law  p mVm  =  p(x)  V  (x)  yields  the  following  instantaneous  pressure  p(x) 

pi*)  -  pm(  1  -  3 xlLy' 

The  pressure  wave  p0(x)  in  the  compression  chamber  is  the  difference  between  this 
instantaneous  pressure  and  the  mean  pressure,  and  is  given  by 

P0(x)  ■  Pi*)  -Pm’P„  «i  -  */3 L)"f  -  l] 


Figure  10  is  plot  of  normalized  gas  spring  force  p0(x)/pm  assuming  isothermal  conditions  (  r  =  1). 
From  this  figure,  it  is  clear  that  the  amplitude  of  pressure  wave  depends  on  the  displacement  x 
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(relative  to  the  length  of  the  compression  chamber).  For  the  prototype,  p0  =  10  psi,  pm  =  30  psi 
yields  p^Pn,  =  0.33.  For  this  data,  this  figure  shows  that  x/L  should  be  0.75  to  achieve  this  value. 
Since  L  -  0.013,  the  diaphragm  should  vibrate  by  x  =  ±  0.010  in.  to  achieve  the  specified 
dynamic  pressure  of  ±  10  psi. 


displacement/length  x/L 


Figure  10.  Gas  Spring  Force  (in  Compression  Chamber)  Fora  Given  Diaphragm 
Displacement. 


Consider  the  Back  Chamber.  The  displacement  x  reduces  the  pressure  from  pm  to  p'(x)  and 
increases  the  volume  from  V0*  to  V(x).  Adiabatic  law  yields  the  instantaneous  back  pressure 
p'(x)  and  associated  back  pressure  wave  p'0  (x) 
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p'C*)  ■  Pm  (1-^/3LTy 


P0(x)  ■  p'(* )  -  pm-  pm[  1  -  xlS X)'Y  -  1] 


This  back  spring  force  is  obviously  a  mirror  image  of  the  gas  spring  force.  The  gas  spring  force  is 
determined  from  gas  pressures  in  the  compression  space  and  back  chamber.  This  force  will 
always  oppose  the  displacement  and  is  directed  towards  the  position  of  static  equilibrium  x  =  0. 

It  is  given  by 


g(*) 


[/>(*)  -  P  '(*)]  A 
(6  xlL)  _  , 


9  -(x/Ly 


Figure  1 1  shows  this  normalized  net  gas  force  g(x)/pmA  as  a  function  of  displacement  under 
isothermal  conditions.  From  this  figure,  it  is  clear  that  the  net  spring  rate  of  the  gas  pressure  is 
nearly  linear.  The  gas  force  can  thus  be  approximated  by  a  gas  spring 

3 p  A 

g(x) - —  x-k  x 

4  L  g 


where 

kg  =gas  spring  rate,  Ib/in. 

pm  =  1 0  psi,  A  =  1 1 .95  in2,  L  =  0.01 3  in.  and  =  1  leads  to  gas  stiffness  of  kg  =  6900  Ib/in. 
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Figure  1 1.  Normalized  Net  Gas  Force  vs.  Displacement 


Since  the  loads  are  axi-symmetric,  the  deflection  varies  with  the  radial  position  only.  The 
deflection  of  the  diaphragm  obviously  depends  on  the  attraction  and  gas  pressures.  They  can  be 
computed  from  formulas  available  from  several  sources  [10-13],  These  formulas  are  based  on 
thin  plate  theory,  i.e. ,  deflections  are  assumed  to  be  smaller  than  the  thickness.  In  our  case  the 
deflection  can  be  larger  than  the  thickness,  so  the  non-linearity  should  be  taken  into  account. 

The  small  deflection  formulas  however  have  the  added  advantage  that  superposition  principle 
can  be  used  to  describe  the  behavior  under  a  variety  of  diaphragm  loading  arrangements.  A 
simple  Fortran  program  was  written  to  obtain  the  deflection  contours,  stresses  and  swept 
volume.  This  program  was  validated  by  hand-calculations  for  several  cases  of  diaphragm 
dimensions  and  loadings. 
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While  testing  the  Phase  I  compressor  model,  we  discovered  a  phenomenon  we  called  "hump 
effect".  This  phenomenon  is  characterized  by  the  inner  region  of  diaphragm  being  pulled  inwards 
by  the  attraction  load  and  the  peripheral  regions  pushed  outward  by  the  gas  loads.  The 
attraction  load  -acting  through  pole  area  -  pulls  the  diaphragm  into  the  compression  chamber. 

On  the  other  hand,  the  compressive  gas  load  -  acting  over  a  larger  diaphragm  area  -  pushes  the 
diaphragm  out  of  the  compression-chamber.  The  deflection  contour  is  the  resultant  of  these  two 
actions.  Because  attraction  pressure  acts  through  the  pole  area  and  the  gas  pressure  acts  over 
the  larger  diaphragm  area,  it  is  possible  that,  for  certain  combinations  of  attraction  and  gas 
pressures,  the  central  portion  of  diaphragm  (facing  the  pole)  will  be  deflected  inwards  while  the 
remaining  portion  might  project  out  of  the  compression  chamber  as  a  "hump".  This  "hump" 
reduces  the  swept  volume.  Obviously  the  attraction  pressure  must  be  sized  to  avoid  this  "hump" 
phenomenon. 

In  order  to  overcome  the  hump  phenomenon,  it  is  necessary  to  accurately  estimate  the  shape  of 
the  deflection  curve  under  widely  varying  loads.  The  diaphragm  deflection  can  be  broken  into  DC 
deflection  -  due  to  DC  loads  and  AC  deflection  -  due  to  AC  loads.  It  is  assumed  that  the  DC 
loads  are  zero,  so  the  DC  deflection  is  zero.  Hence  we  focus  on  AC  deflections  due  to  AC  load 
only.  These  loads  are  given  by  equation  (5).  The  superposition  principle  is  used  to  calculate  the 
net  deflection  under  these  variety  of  loads.  In  general,  the  central  deflection  is  estimated  from 

k  x  -  f 

where 

k  =  stiffness  of  diaphragm  under  load  f 
x  =  deflection  of  diaphragm 
f  =  AC  load  on  the  diaphragm 

The  stiffness  of  diaphragm  can  be  estimated  from  known  formulas  such  as  those  given  in 
Roark. 

Let  us  assume  that  the  compression  and  back  chambers  are  identical.  The  central  deflection 
due  to  attraction  loads,  which  contain  the  AC  current  component  and  deflection  component,  is 
computed  from 
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k  X 
a 


-  2k[gc  -  2k  x  i.e., 
(kd  -  2k)x  =  2  k/gc 


where 

ka  =  stiffness  of  diaphragm  under  attraction  load 

From  this  it  is  clear  that  the  net  stiffness  is  the  diaphragm  stiffness  minus  twice  the  negative 
stiffness.  Since  net  stiffness  is  softer  than  the  diaphragm,  smaller  current  will  be  required  by  the 
electromagnet. 

The  central  deflection  due  to  the  gas  loads  are  similarly  calculated  using  gas  springs.  Note  that 
the  effect  of  gas  springs  is  to  increase  the  net  stiffness  whereas  the  effect  of  attraction  stiffness 
is  to  reduce  it.  Accurate  estimation  of  stiffness  of  these  springs  is  therefore  vital  to  determining 
the  central  deflection  of  the  diaphragm. 


The  dynamic  or  AC  stresses  in  the  diaphragm  are  of  significant  interest  as  they  strongly 
influence  the  fatigue  life  of  the  machine.  The  phase  difference  between  the  pressure  wave  and 
deflection  also  influences  the  stress  amplitudes.  In  Phase  I  however,  to  simplify  the  analysis, 
we  assume  that  the  attraction  pressures  are  180  degrees  out  of  phase  with  the  gas  pressure 
wave.  During  one  complete  cycle,  the  diaphragm  will  fully  deflect  into  the  back  cavity  plate, 
reverse  itself  and  fully  deflect  into  the  compression  cavity  plate.  The  stresses  in  the  diaphragm 
are  the  membrane  radial  stress,  membrane  tangential  stress  and  the  bending  tangential  stress. 
These  stresses  are  then  combined  algebraically  to  yield  maximum  stress  which  governs  the  life 
of  the  diaphragm.  The  combined  stress  is  maximum  at  the  inflection  point  and  at  the  center  of 
the  cavity.  The  inflection  point  is  the  point  where  the  two  radii  are  tangent  to  each  other.  For 
small  deflections  of  thin  plate  loaded  by  pressure  q  spread  from  radius  r  =  r0  to  r  =  a,  the  radial 
and  tangential  stresses  can  be  computed  from  the  formulas 
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The  deflections  are  usually  larger  than  half  the  thickness  of  the  diaphragm.  In  such  cases,  the 
middle  surface  becomes  appreciably  strained  and  the  stress  in  it,  called  the  membrane  stress, 
cannot  be  ignored.  This  stress  enables  the  diaphragm  to  carry  a  part  of  the  load  in  direct 
tension. 

Under  these  conditions,  the  maximum  stress  can  be  estimated  from  (see  Roark,  p.  433) 

(4.267  tx  ♦  0.476  x2) 


where 

a  =  maximum  stress  in  the  diaphragm 
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i  =  thickness  of  diaphragm 
x  =  maximum  deflection  of  the  diaphragm 
E  -  Young's  Modulus 


4.. 5  Diaphragm  material 

The  diaphragm  material  must  be  carefully  selected  to  meet  the  long  life  requirement.  The  most 
popular  materials  are  Ti-6Ai-4V  and  stainless  steels.  Major  considerations  in  selecting  the 
diaphragm  material  are: 

o  it  must  operate  without  failure  over  1010  fatigue  cycles 

o  it  should  be  capable  of  developing  large  deflections  without  failure 

o  it  must  not  react  with  Helium 
o  its  outgassing  must  be  very  low 

o  it  should  be  available  in  the  thin  sheets  0.010  -  0.080  in. 

o  it  must  be  light  weight  (high  strength  to  wt.  ratio). 

In  order  to  have  long  life,  ideally,  a  diaphragm  should  have 

o  high  yield  strength  to  yield  the  long  life  requirement 

o  low  Young’s  Modulus  to  meet  the  large  swept  volume/pressure  wave  requirement 


A  high  yield  strength  ensures  that  the  diaphragm  will  carry  large  stresses  without  failure.  A  low 
Young's  modulus  on  the  other  hand  makes  the  diaphragm  flexible  enough  to  create  a  large 
swept  volume.  A  high  yield  strength/Young's  modulus  will  therefore  signify  longer  life  while 
delivering  larger  swept  volume.  The  ratio  of  yield  strength  to  Young's  modulus  is  called  the  life 
factor  oE. 

Yield  Strength  o 

aE  -  - — 

Young  's  Modulus  E 
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The  central  deflection  is  proportional  to  the  ratio  of  maximum  stress  and  Youngs  modulus. 

Since  higher  deflection  obviously  means  higher  swept  volume  and  hence  higher  pressure  wave, 
to  maximize  the  pressure  wave,  high  life  factor  is  desirable.  Thus,  the  most  appropriate  material 
for  the  diaphragm  is  one  which  has  the  largest  life  factor. 

Recognizing  the  importance  of  the  life  factor,  we  have  collected  data  on  this  value  for  various 
materials.  These  materials  may  be  broadly  classified  into  metallic  and  non-metallic.  We  note  that 
nonmetallic  materials  can  be  as  strong  as  metallic  materials.  But  we  eliminated  them  for  space 
cryocooler  application  as  they  have  significantly  higher  outgassing  than  metallic  materials.  Table 
1  below  indicates  the  life  factor  of  typical  metallic  and  non  metallic  materials. 


Table  1 .  Long  Life  Factor  oE  for  Membrane  Materials  (at  300  K) 


N 

0. 

Material 

Density 

lb/in3 

Yield 

Strength 

aM  ksi 

Endur¬ 
ance  Limit 

0.„„  ksi 

Young's 

Modulus 

E  10'psi 

Long  life 

Factor 

0F  =  O^E 

Ref. 

i. 

Ti-6AI-4V 

Annealed 

0.160 

1207 

80* 

14-19 

8.6 

5-6 

2 

Chromium  steel 

7C27M02 

0.3 

210 

210/2  = 

100 

30 

7 

7 

2. 

Maraging  Steel 

18Ni(300) 

0.3 

280 

280/2  = 

140b 

27 

10.4 

3. 

Be-Cu  172, 1-2 

Hard 

0.3 

80* 

18 

4.4 

6 

4. 

Stainless  Steel 

321  Annealed 

0.29 

42* 

28 

1.5 

6 

5. 

Aluminum  6061- 

T6 

0.097 

30* 

14c 

10 

6. 

Polyamid  Nylon  6 

Comturf  613 

.048 

272 

1 35b 

1.0 

272 

1 

Polyster  Glass 

Fibre  Reinf. 

Hyside  1 1  YRr. 

0.068 

240 

8. 

Polyster  Woven 

Glass  Cloth 

80 

406 

2.5 

32 
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a)  At  77K  (see  Ref.  11) 

b)  Considered  as  half  of  flexural  yield  strength 

c)  see  ref  1 1 

7)  Guaranteed  min.  per  catalog  R  President  Titanium,  Hanson,  MA  02341. 


From  this  table  it  is  clear  that  Titanium  alloy  Ti-6AI-4V  and  maraging  steel  are  ideal  metallic 
materials.  The  beryllium  copper  has  a  lower  life  factor  than  the  two  and  hence  is  not  preferred. 
Since  we  need  a  magnetic  material  for  the  diaphragm  (in  order  to  attract  the  electromagnet),  we 
believe  that  maraging  steel  is  an  ideal  candidate.  However,  this  steel  is  not  available  in  sheet 
form  so  we  prefer  the  next  best,  viz.  chromium  stainless  steels  as  possible  candidate. 


5.2.6  Permissible  Working  Stress 


A  membrane  vibrating  at  60  Hz  for  10  years  undergoes  2  x  1010  cycles  of  reversal  of  stress. 
Fatigue  data  at  these  large  number  of  cycles  is  not  readily  available.  Hence  we  rely  on  the 

Yield  Strength  a 

N  ,  - - - 2" - 

Max .  Working  Stress  o 

°  max 


"Fatigue  Safety  Margin  N"  concept.  It  is  defined  above.  The  yield  strength  is  taken  at  0.2  % 
amplitude.  The  maximum  working  stress  is  the  amplitude  of  cyclic  stress  that  will  ensure  fail¬ 
safe  operation  up  to  2  x  1010  cycles.  As  rule  of  thumb,  we  take  N  =  10  or  greater  to  ensure  life 
of  ten  years  as  per  the  rationale  given  below. 

Wong  et  al  [3]  indicates  that  for  a  high  quality  stainless  steel,  spring  material,  the  peak 
stress  should  be  limited  to  one-half  of  the  endurance  limit,  or  415  MPa  (60,000  psi)  in 
order  to  achieve  the  10  year  life  requirement.  Since  endurance  limit  is  usually  one-half  of 
yield  strength,  the  Fatigue  Safety  Margin  according  to  him  is  N  =  4. 

In  contrast,  Stacy  [11]  indicates  that  the  membrane  stresses  should  be  less  than  90  MPa 
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(13,000  psi)  to  achieve  long  life.  In  fact,  Stacy  says  that  100  MPa  peak  stress  limit  falls 
well  below  99.999%  fatigue  probability  at  1010  cycles  of  recurring  stress.  The  yield  stress 
is  120  ksi  for  the  Ti-6AI-4V  material  he  used.  The  Fatigue  Safety  Margin  is  N=  9. 

Thus,  there  is  a  significant  difference  in  the  Fatigue  Safety  Margins  specified  by  the  two  experts. 
Therefore,  as  a  matter  of  abundant  precaution,  we  choose  N  =  10.  This  yields  a  maximum 
working  stress  of  10,000  psi  for  a  Ti-6  AI-4V  material.  Our  decision  to  limit  stresses  to  10000  psi 
will  exceed  both  above  safety  margins  and  hence  will  have  long  life. 

A  classic  fatigue  curve,  such  as  that  shown  in  Figure  12  shows  that  as  the  stress  level  is 
increased,  the  number  of  cycles  to  failure  decreases.  Conversely,  as  the  stress  level  is 
decreased,  the  number  of  cycles  to  failure  increases.  Typically,  a  stress  level  of  100  ksi  results 
in  a  life  of  100  to  1000;  a  stress  level  of  75  ksi  will  yield  a  life  of  104  to  105  cycles.  Finally  by 
dropping  the  stress  level  to  75  ksi,  the  life  goes  up  to  10®.  Reducing  it  further  to  about  20  ksi  will 
increase  the  life  further  to  1010. 

Consider  a  strain  hardened,  cold  rolled  301  SS  material  whose  properties  are  far  in  excess  of 
annealed  material.  This  material  however  is  non-magnetic  and  hence  could  not  be  selected  for 
our  compressor.  The  yield  strength  of  this  material  is  180,000  psi.  A  fatigue  life  drawing  such  as 
figure  indicates  that  a  design  stress  below  40  ksi  can  reach  life  of  1010  cycles. 
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THESE  CURVES  ARE  BASED  ON  DIAPHRAGM  MATERIAL 
FOLLOWING  MECHANICAL  PROPERTIES: 


ULTIMATE  TENSILE  STRENGTH 
YIELD  STRENGTH 
ELONGAT I  ON 
REDUCTION  OF  AREA 
HARDNESS 


-  190.  000  PS  I, 

-  180,  000  PSI, 

-  15 %  MIN. 

-  35  TO  44% 

-  38  TO  45  RC 


MIN. 

MIN. 


5.2.7  Fatigue  Analysis 

Since  the  membrane  fails  by  flexure,  the  flexural  yield  strength  can  be  taken  as  a  rough  measure 
of  the  fatigue  strength.  However,  this  data  is  not  readily  available  for  many  materials,  so  yield 
strength  is  taken  as  a  primary  measure  of  fatigue  strength. 

Since  diaphragm  stresses  are  cyclic,  the  cyclic  stress  must  be  addressed  during  the  design 
stage  itself.  The  fatigue  data  can  broadly  be  classified  into  low  cycle  data,  which  is  determined 
by  strain  cycling  to  a  level  that  is  beyond  the  yield  strength.  A  common  example  is  tensile  test  to 
determine  the  yield  strength  elongation  etc.  A  high  cycle  fatigue  data  is  accumulated  while  the 
material  is  being  stressed  in  the  elastic  strain  region. 


In  addition  to  the  stress,  factors  that  influence  the  fatigue  life  are: 
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o  Surface  finish  -  The  surface  finish  of  the  diaphragm  must  be  highly  polished  and  smooth  to  be 
free  of  nicks,  scratches  and  metal  inclusions  in  order  to  attain  high  cycles.  Surface  scratches 
and  nicks  produce  stress  concentrations  of  varying  magnitude  ranging  from  1  to  3  and  the  stress 
level  will  increase  rapidly. 

o  Tensile  Strength  -  Higher  tensile  strength  usually  yields  higher  fatigue  life.  However,  when 
tensile  strength  is  increased,  the  material  becomes  very  notch  sensitive,  and  a  point  is  reached 
when  additional  strength  compromises  toughness. 

o  Thickness  -  The  stresses  are  inversely  proportional  to  the  thickness  to  produce  a  given 
pressure  wave.  When  the  thickness  is  less  than  0.010  in.,  the  material  is  difficult  to  handle  and 
cannot  stand  negative  stiffness.  Typical  diaphragms  therefore  range  from  0.015  to  0.020  in 
thickness.  Other  things  being  equal,  a  thinner  material  will  fail  sooner  than  thicker  material  due 
to  crack  propagation. 

o  Cavity  shape  -  The  shape  of  the  cavity  exerts  the  greatest  influence  on  the  stress  level.  The 
cavity  must  be  precisely  machined,  preferably  using  a  numerically  controlled  lathe  with  the  cavity 
automatically  generated  by  a  program.  Once  the  cavity  is  cut,  it  must  be  inspected  using  a  cavity 
plotter  to  verify  the  contour  within  0.001  in.  Maintenance  in  the  field  can  also  be  an  important 
factor.  Excessive  polishing  or  lapping  the  head  could  unfavorably  change  the  stress  levels. 

o  Torque  -  The  head  bolts  must  be  torqued  in  a  cross  pattern,  incremental  and  to  specific  levels. 
If  the  edges  of  the  diaphragm  are  not  uniformly  clamped  and  restrained,  the  stresses  in  it  change 
significantly. 

o  Cleanliness  -  Dirt  in  any  form  will  change  the  cavity  shape  and  cause  stress  concentrations. 
The  smallest  of  the  particles  that  enter  into  the  cavity  will  usually  work  their  way  to  the  outer  edge 
and  build  up  on  the  surface  adjacent  to  the  clamped  seal  area.  Larger  particles  will  cause  stress 
concentrations  and  can  change  stresses  by  a  factor  of  1.2  to  3. 

o  Fretting  Corrosion  -  This  fretting  occurs  when  there  is  low  frequency  relative  motion  between 
two  elements  under  load.  Metal  is  removed  from  one  element  and  this  sets  up  stress 
concentration,  usually  in  the  diaphragm.  To  reduce  the  fretting,  the  contact  surfaces  must  be 
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relieved.  In  addition,  high  temperatures  will  affect  the  diaphragm  stresses  significantly. 

301  Steel.  A  strain  hardened,  cold  rolled  301  SS  material  has  properties  that  are  far  better  than 
that  of  annealed  material.  The  yield  strength  of  this  material  is  180,000  psi.  A  fatigue  life 
drawing  such  as  figure  12  indicates  that  a  design  stress  below  40  ksi  can  reach  life  of  1010 
cycles.  This  material  is  non-magnetic  and  hence  could  not  be  selected  for  our  compressor. 

Sandvik  7C27M02  steel.  This  steel  sheet  has  a  tensile  strength  of  261,000  psi.  It  has  a  fatigue 
limit  of  125,000  psi.  Limiting  the  stress  to  10,000  psi  will  therefore  greatly  increase  the 
confidence  that  it  will  work  without  failure. 


4 . 6  Dynamic  SViodiel  of  the  Diaphragm 

The  diaphragm  can  be  modeled  in  a  number  of  ways.  In  the  following  f  denotes  the  net  attraction 
force  and  g  denotes  the  net  gas  force.  The  first  model  is  the  static  model  kx  =  f0-g0 ,  [used  by 
Creare,  Stacy  and  Shimko  2]  which  ignores  the  dynamic  effects  and  assumes  that  the 
diaphragm  deflects  statically.  The  next  refinement  is  the  pressure  wave  model  m  x  +  cx  +kx  = 
f(t)-g(t)  which  assumes  that  attraction  and  gas  pressure  do  not  change  with  displacement.  This 
model  was  used  by  Marquardt,  Radebaugh,  Kittel  [9],  A  further  refinement  is  the  gas  spring 
model  mx  +  cx  +  kx  +  g(x)  =  f(x)  which  recognizes  that  the  attraction  and  gas  forces  behave  as 
non-linear  spring.  This  model  recognizes  that  gas  pressures  change  with  the  displacement. 

This  model  was  used  by  Sussholz  [10]  for  pressure  sensors  and  Jonge  [11]  and  Boyle  et  al  [12] 
for  piston-driven  compressors.  We  will  investigate  this  mode!  in  Phase  II. 

Note  that  in  Phase  I  we  deal  with  an  unloaded  compressor.  The  model  is  not  connected  to  a 
displacer.  As  such  the  pressure  wave  that  is  generated  will  be  confined  to  the  compression 
chamber.  This  pressure  wave  is  not  transmitted  into  a  displacer.  In  this  sense  the  Phase  I 
compressor  does  not  realistically  depict  its  performance.  In  the  real  world,  the  compressor  is 
always  connected  to  a  displacer  using  a  small  tube.  The  volume  of  the  displacer  chamber  adds 
to  that  of  the  compression  chamber  and  this  effects  the  strength  of  the  pressure  wave.  In  Phase 
II  we  will  develop  the  more  realistic  model  of  the  loaded  compressor  which  includes  the  effect  of 
volume  of  displacer. 
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We  make  following  assumptions  in  developing  the  model:  (1)  The  pressure  within  the 
compression  chamber  follows  the  law  of  adiabatic  compression  [13],  PV  =  constant,  =1.4  (2) 
The  dissipative  forces  such  as  heat  loss,  friction  and  leakage  are  small,  and  can  be  modeled  by 
a  linear  viscous  damper.  (3)  The  internal  pressures  are  uniform  throughout,  and  diffraction 
effects  at  the  mouth  of  the  outlet  are  negligible.  (4)  The  continuous  diaphragm  can  be  modeled 
by  an  equivalent  discrete  mass-spring  system. 

From  the  free  body  diagram  shown  in  Figure  13,  the  displacement  x  is  governed  by 

mx  +  cx  +  kjc=f-g(x) 


where 

m  =  effective  mass  of  the  diaphragm,  Ib/g 
c  =  effective  damping  coefficient  of  the  system,  Ib-sec/in 
kd  =  effective  stiffness  of  the  diaphragm,  Ib/in 

x(t)  =  Xc  sin  ( oot-cj))  =  instantaneous  displacement  of  diaphragm's  center,  in 
g[x(t)]  =  gjx]  sin  (go  t-ip )  =  restoring  gas  force  on  the  diaphragm 
f[x(t)]  =  f0[x]  sin  go  t  =  applied  attraction  force  on  the  diaphragm 

Initially,  there  is  no  attraction  force;  compression  and  back  chamber  pressures  are  equal  and 
opposite.  So  the  diaphragm  is  in  equilibrium  at  x  =  0. 


z  a , 


Figure  13.  Free  Body  Diagram  of  the  Diaphragm. 
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4 . 6 .  *1  Estimation  of  Attraction  Force 


The  attraction  pressure  pa  moves  the  diaphragm  to  the  right,  which  increases  the  fluid  pressure 
in  the  right  captive  chamber  from  pm  to  pm  +  p0sin(a)t-ijj)  .  It  also  reduces  pressure  in  the  left 
chamber  from  pm  to  pm  -  p,  sin(wt-  ip) .  The  dynamic  pressures  p0  and  p,  push  the  diaphragm  to 
the  left.  In  other  words,  the  attraction  pressure  is  always  opposed  by  the  dynamic  pressures  it 
creates.  Let  Xc  denote  the  central  deflection  of  the  diaphragm  under  the  combined  action  of  pa, 
p0  and  p,.  The  swept  volume  Vs  during  peak-to-peak  displacement  is  the  twice  conical  volume 
contained  by  base  diameter  2a  and  the  central  deflection  x,.,  and  is  given  by 

F  .  -k  a  V 
*  3 

The  dynamic  pressure  amplitude  is  related  to  this  swept  volume  by: 

£*  Vsji_ 

Pm  "  & 

m 


Substituting  (2)  in  (3)  shows  that  the  central  deflection  can  be  computed  from  specified  pressure 
wave  amplitude  by 

Pr 

x  0  - 

C  «  a2/ 3 


44 


5.  COMPRESSION  CHAMBER 


5.1  Design  Requirements 

The  compression  chamber  provides  a  dirt-free  space  for  creating  a  specified  pressure  wave. 

The  volume  between  a  vibrating  diaphragm  and  the  contour  of  cylinder  defines  this  compression 
space.  The  design  objective  is  to  size  this  volume  from  the  following  data: 

pm  =  mean  charge  pressure 

r  =  pressure  ratio 

a  =  ratio  of  deflection  volume  and  swept  volume 
AV  =  swept  volume 

Vp  =  payload  volume 

y  =  specific  heat  ratio 

In  addition,  we  assume  that  the  diameter  D  of  the  compression  chamber  is  specified  a  priori. 

The  mean  captive  volume  is  the  initial  volume  into  which  the  working  fluid  is  charged.  It 
comprises  of  the  pay  load  volume  (displacer  volume)  and  the  compression  chamber  volume. 
The  net  volume  .  It  is  broken  into  a  deflection  volume  and  a  dead  volume.  The  deflection 
volume  is  the  volume  between  the  undeflected  and  deflection  positions  of  the  diaphragm.  The 
dead  volume  is  the  balance  of  mean  captive  volume  that  is  left  after  the  diaphragm  fully  deflects. 
The  volume  parameters  to  be  determined  are  hence: 

V  =  mean  captive  volume  (=  Vp+ Vw) 

Vw  =  compression  chamber  volume  (=  VD  +  Vd) 

VD  =  deflection  volume  (  =  aAV) 

Vd  =  dead  volume  (=  Vw  -  VD) 

The  adiabatic  gas  law  PV  =  const,  will  be  used  to  size  the  compression  chamber.  This  law 
defines  the  behavior  of  an  ideal  gas  when  its  volume  changes  with  no  heat  entering  or  leaving. 
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This  condition  is  usually  satisfied  during  the  compression  phase  when  pressure  changes  from  Pm 
to  Pmax  and  volume  changes  from  v  to  Vmin .  [  The  diaphragm  movement  gives  rise  to  gas  flow 
causing  the  pressure  difference  and  the  displacer  position  determines  the  amount  of  gas  in 
compression  and  expansion  spaces;  these  effects  are  ignored  in  our  elementary  model] 

P  max  ^ m  If  ,  . 

-  o  -  e  -  (a) 


Noting  that  minimum  and  mean  captive  volumes  are  given  by: 


v_ 


v  . 


f 


V  -  a  AV 

W 


(b) 


Using  (b)  in  (a)  yields  following  expression  for  compression  chamber  volume: 


-a  —  AV  -  V  where  R  =  ( — 

R  -  1  p  r  *  l 


\i/y 


(1) 


This  compression  chamber  volume  is  used  to  compute  the  volume  parameters  using: 

(2) 


V  -  V  *  V 

m  p  \? 

VD  =  a  A  V 

VD 


Axial  Length  of  Compression  Chamber 


The  axial  length  of  a  cylindrical  compression  chamber  is  computed  from: 


-  D  2L  -  F  i.e,  L 


4F. 
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where  L  =  axial  length  of  a  noncontoured  compression  chamber. 

If  this  L  is  found  to  be  less  than  maximum  central  deflection,  then  a  contoured  compression 
chamber  will  be  required.  Using  a  cone  approximation  formula,  the  axial  length  of  a  contoured 
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compression  chamber  will  be  found  from: 


iiiii.f,  u, 

3  4 


12  F 
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5.2  Back  chamber 


The  back  chamber's  primary  function  is  to  balance  the  diaphragm  against  the  charge  pressure. 
This  charge  pressure  pushes  the  diaphragm  outward  and  creates  static  stresses  in  it.  There  are 
two  methods  to  counter  this  action  of  charge  pressure: 

(a)  Sealed  Back  Chamber  Approach.  (Figure  14a).  This  approach  uses  two  attraction 
electromagnets,  one  is  the  compression  chamber  and  another  is  the  back  chamber  (i.e.,  it 
requires  electromagnets  on  both  sides  of  the  diaphragm).  The  back  chamber  is  charged  with 
gas  at  mean  pressure  and  is  sealed.  The  diaphragm  thus  remains  at  the  neutral  position  under 
the  action  of  two  equal  and  opposite  gas  pressures  as  shown  in  this  Figure.  However,  this 
approach  has  the  disadvantages  that  it  creates  a  pressure  wave  in  the  back  chamber  (that  is  not 
used  in  the  cryocooler).  This  unused  pressure  wave  opposes  the  displacement  of  the 
diaphragm.  The  attraction  magnet  should  therefore  overcome  this  unused  pressure  in  addition 
to  the  useful  pressure  wave  in  the  compression  chamber.  So  this  approach  increases  external 
force  required  without  increasing  the  output  compression.  As  such,  this  approach  is  not  a 
preferred  method. 

(b)  Open  Back  Chamber  Approach  (Figure  14b).  This  uses  a  single  electromagnet  on  the 
compression-chamber.  It  does  not  use  any  back  chamber  (i.e,  one  side  of  the  diaphragm  is  left 
open  to  atmospheric  pressure).  It  uses  a  permanent  magnet  bias  to  neutralize  the  DC  stress 
produced  by  the  charge  pressure.  This  permanent  magnet  can  be  positioned  in  the  cylindrical 
face  of  the  compression  chamber  as  shown.  As  a  result  the  diaphragm  initially  deflects  inwards 
into  the  TDC.  The  compression  chamber  is  then  charged  with  pressurized  gas  to  bring  it  back  to 
the  mean  position.  The  magnet  is  sized  so  that  its  attraction  force  balances  the  charge  pressure 
and  positions  the  diaphragm  at  its  neutral  position. 

The  advantage  of  this  approach  is  that,  unlike  the  sealed  back  chamber  approach,  it  does  not 
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create  an  additional,  but  unused  pressure  wave  in  the  back  chamber.  As  a  result,  this  approach 
requires  less  attraction  force  than  that  needed  by  the  sealed  back  chamber  approach.  The  use 
of  a  smaller  attraction  force  obviously  saves  significant  input  power,  and  hence  this  approach  is 
the  preferred  approach  .  This  approach  will  be  investigated  in  the  Phase  II  of  the  project. 

To  qualify  the  advantage  of  the  open  back  chamber  approach  we  considered  the  problem  of 
generating  a  25  ±  8.33  psi  pressure  wave.  We  used  a  3.9  in  diameter  0.015  thick  steel 
diaphragm  to  generate  this  pressure  wave.  The  sealed  back  chamber  approach  produced  an 
additional  ±  8.33  psi  pressure  wave.  Our  computations  showed  that  a  ±  30  psi  attraction 
pressure  (at  1.0  pole  diameter)  will  be  required  to  overcome  the  net  ±16.67  psi  pressure  wave  on 
the  diaphragm  if  the  sealed  back  chamber  approach  is  used. 

If  the  open  back  chamber  approach  is  used,  the  attraction  pressure  needs  to  overcome  ±  8.3  psi 
of  pressure  wave  in  the  compression  chamber  only.  As  a  result  our  analysis  indicates  that  the 
attraction  pressure  can  be  reduced  from  ±30  psi  to  ±  15  psi.  As  a  result  we  save  half  the  power 
needed  to  produce  the  desired  pressure  wave. 
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(b)  PermanentMagnet  Approach 


Figure  14.  Two  Design  Approaches  to  Back  Chamber  Design. 


This  report  presented  the  investigations  carried  out  by  Precision  Magnetic  Bearing 
Systems  under  an  SB1R  project  sponsored  by  Ballistic  Missiles  Defense  Organization. 
The  report  addressed  the  issue  of  feasibility  of  a  novel  electromagnetic  diaphragm 
compressor  to  drive  cryocoolers.  Two  prototype  compressor  units  were  built  to 
establish  the  feasibility.  The  first  prototype  used  permanent  magnets  to  bias  the 
vibrating  diaphragm  while  the  second  prototype  used  electromagnets.  The  diaphragm 
used  in  the  first  prototype  was  found  to  be  too  stiff  and  hence  electromagnet  was  used 
in  the  second  prototype  together  with  a  thinner  diaphragm.  Testing  of  the  second 
prototype  indicated  that  it  is  possible  to  develop  pressure  ratios  of  the  order  of  1.67  at 
frequencies  of  up  to  60  Hz.  Since  cryocoolers  require  pressure  ratios  of  about  2  at 
frequencies  of  about  40  Hz,  it  is  concluded  that  it  is  feasible  to  develop  an  engineering 
model  of  the  compressor  to  drive  cryocoolers. 
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